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00 ' ABSTRACT 

oo : 

f"**) . Context. We present work on Ha spectral line characteristics in PHOENIX stellar model atmospheres and their comparison to microlensing 
"^" ' observations. 

<^^> Aims. We examine in detail the Ha equivalent width (EW) and the line shape characteristics for effective temperatures of 4500K< T e g < 5600K 
>0 ■ where Ho- is a strong spectral feature. 

Methods. We find that Ho- EW in models calculated under the assumption of local thermodynamic equilibrium (LTE) is up to 15% smaller 
ff , than in models without this assumption, non-LTE models (NLTE) and that line shapes vary significantly for the two model types. A comparison 
O i 1 with available high quality microlensing data, capable of tracing Ho- absorption across the face of one G5III giant, shows that the LTE model 
that fits the EW best is about 100K hotter than and the best-fitting NLTE model has a similar T eff as predicted by the spectral type analysis of 
the observed star but agree within the uncertainties of the observationally derived temperature. 

Results. Neither LTE nor NLTE models fit the line shape well. We suspect unmodelled chromospheric emission. Line shape diagnostics 
frt suggest lower gravities than derived for the star and are unacceptable low in the case of the LTE models. We show that EW alone is insufficient 
• • • for comparison to stellar model atmospheres, but combined with a new shape parameter we define is promising. In stellar parameter ranges 
. ^- w h ere the Her line is strong, a NLTE approach of modeling stellar atmospheres is not only beneficial but mandatory. 
K> : Conclusions. 

5^ , Key words, stars: atmospheres - modeling - Gravitational lensing: microlensing 

1 . Introduction biggest advance in one-dimensional (1-D) atmospheres model- 

ing has been to introduce consistent treatment of Non-Local 
Mostof our basic knowledge about the interior physics of stars, r™ ■ • c .,., ■ , XIITn 1 A , I 

& F J Thermodynamic Equilibrium (NLTE) processes (Anderson, 

even of the Sun, is based on our understanding of stellar at- linorlltj TTTTcd "llinnrlv u , .• , . . 

6 1989; Hauschildt & Baron, 1999). However observational tests 

mospheres. The modeling of such atmospheres allows us to fMITr a . . , ,.„, ... . ... , , , , 

r or f NLTE effects have been difficult to establish, because study- 
derive the star's fundamental characteristics such as temper- . , c . .. . ... c , 

r ing atmospheres of stars more distant than our Sun has proven 
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ature, gravity, chemical composition and age. However many ,. ffi , , ,. 
outstanding questions still need to be answered about stellar 
atmospheres. The "solar model problem," , i.e., the disagree 



Various methods have been tried to resolve stars, 



ment between model predictions for the neon-oxygen abun- ie ' measuring their sizes and resolving surface struc- 
dance ratio and measurements by helioseismology, is hotly de- ture ' Direct observations have been restricted to the Sun 
bated (Drake & Testa, 2005; Asplund et al., 2003). Derivation ^' Blackwe11 - Lvnas-Gray. & Smith, 1995) or, when us- 
of elemental abundances is based on the accuracy of stellar in 8 the Hubble Space Telesc opefflST), to nearby super- 
model atmospheres and as such make and shake our knowl- « iants (e «' Gilliland & DuDree ' 1996) ' Do PP ler ima S in S ^ 
. , , . . ,, , . , , . . .. , ■ Olah, Jurcsik, & Strassmeier, 2003) produces observed center- 
edge not only about stellar, but also about galactic and cosmic ■ ■ ■ •■ " F 

evolution (As D lund,|Sl>. to - limb variations of the integrated flux that are in good 

rj., ,.rr- ,. . t ii u j agreement with values derived from current LTE model at- 

The main difficulty in constraining stellar atmosphere mod- ° I — ■ 1 

, . , , . ., . . . , . , mospheres. Whereas Interferometry (e.g. Aufdenberaet al., 

els with observations is that most techniques only measure i 1 J ° ' ■ " 

...-,■ . . a j- i • i j c . t-u 2005) show less limb darkening than predicted by stellar at- 

the disk-integrated flux or disk-integrated spectra of stars. The ■ v o ±- j 

mospheres models. Due to the uneven surface of the moon, 

Send offprint requests to: C. Thurl surface measurements of nearby giants via lunar occulta- 
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tions (IRichichi & LisiL Il990t) have not yet yielded sufficient 
precision to d istinguish betwe en various model atmosp heres. 
lO'Donoehue et alJ <2003l> and IPopper & Guinanl dl998l) have 
used eclipsing binaries to study surface structure. In such sys- 
tems both stars are often in a similar evolutionary state, which 
complicates tremendously the decomposition of the single star 
spectra. Furthermore, if one star is evolved and has (almost) 
filled its Roche Lobe so that mass transfer occurs, then neither 
star is a good probe for single-star atmospheres. However all of 
these methods are restricted to nearby stars. 

Microlensing was proposed to offer an elegant solu- 
tion to the problem of observing distant stellar atmospheres 
(Schneider, Ehlers, & Falco, 1992). During a microlensing 
event a differential magnification pattern across the face 
of the source star enables the observer to resolve its sur- 
face (see ^5] for more detai l). Several authors hav e stud- 
ied the effect theoretically. iHevrovskv & Sasselovl (J2000) 
studied how microlensing ca n be used to detect ste llar 
spots on the surfaces of stars. Hevrovskv. Sasselov. & Loeb 
(2000) demonstrated the effects of microlensing on syn- 
thesized optical spectra of red giant model atmospheres. 
Hevrovskv (2003) discussed the difficu lty of microlensing^ 
light curve inversion methods. Brvce, Hendrv, & Valls-Gabaud 
J2002l) and iHendrv. Brvce. & Valls-Gabaudl ( J2002h examined 
the detectability of star spots via microlensing and their 
possible effects on mu l ticolor microlensing light curves. 
Brv ce. Ignace. & Hendrvl (120031) examined the change in line 
profiles due to bulk motion in circumstellar envelopes during 
microlensing fold caustic crossin g even ts. 

An early attempt by iLennon et ail ill 9961 1 19971) used the 
magnification effect of the microlensing event 96-BLG-3 for 
spectral studies of stellar atmospheres. They further suggested 
an effect of center to limb variations on the Ha line profile 
during such an event. A number of well-studied microlensing 
events have resulted in measurements that resolve the source 
star; we discuss one of these in more detail in 35.11 The names 
and major characteristics of these events are listed in Table 
n with the respective publications. The major findings were 
promising comparisons of observationally-deduced limb dark- 
ening coefficients with model atmospheres and variations of 
equivalent width (EW) of the Ha line. One of these events, 
EROS BLG-2000-5, is a cusp-crossing event, which shows 
that despite the increased complexity in modeling not only 
fold caustic crossings can be used for thes e purposes. For the 
mos t recent event, OGLE-2002-BUL-69 (Cas san et all 12004: 

1 I — — — i — — 

Kubas et al., 2005), a well-sampled light curve and high res- 
olution spectra are available at crucial times during the event, 
making it a good target with which to compare high resolution 
model atmospheres. Four events of similar quality have been 
observed in 2004 and another four in 2005 and many more 
are to be expected in the years to come, making microlensing a 
viable future for testing stellar model atmospheres. 

This is the first in series of papers aimed at developing 
methods to test stellar atmosphere models by confronting them 
with observations that resolve stellar surfaces photometrically 
or spectroscopically. Based on the most current theoretical 



Event Name 



Publication 



results 



MACHO Alert 95-30 
MACHO 97-BLG-28 
MACHO 97-BLG-41 
OGLE-1999-BUL-23 
MACHO 98-SMC-l 
EROS BLG-2000-5 



MOA 2002 -BLG-23 
OGLE-2002-BUL-069 



Alcocketal. (1997) 
Albrow et al. (1999a) 
Albrow et al. (2000) 
Albrow et al. (2001a) 
Albrow et al. Q_999b) 
Albrow et al. (2001b) 
Castro etal. (2001) 
Fields et al. (2003) 
Abe et al. (2003) 
Cassan et al. (2004) 



Ho-EW 
limb darkening 
limb darkening 
limb darkening 
limb darkening 
limb darkening 
HaEW 

limb darkening 
HaEW 
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Table 1. Observed microlensing events with characteristics and 
results. In the first column are listed microlensing events dis- 
cussed in greater detail by the authors in the second column. 
The last column lists the main findings for the events, i.e. 
Ha equivalent width variations and limb darkening coefficients 
were measured in the respective events. 



models, we describe what observations can be made in the fu- 
ture in order to distinguish between various stellar models and 
learn the most from comparisons to observations. 

We lead into this work by describing, in |3 our methods 
and the PHOENIX stellar models. In ij3]we explain the reasons 
for and analysis techniques of studying the Ha line. In |4] we 
introduce the spectral parameter range studied and show that 
spectral differences between LTE and NLTE cannot be ignored 
within the Ha line for stars with 4500/: < r eff < 5600/:. How 
microlensing can be used as a tool to resolve the surfaces of 
stars is described in |5] where we also present results of our 
comparison of model atmospheric data to microlensing Ha EW 
line data for the event OGLE-BULGE-2002-069. 

In |6] we illustrate how to most appropriately choose the 
best-fitting model out of a set of models for a given parameter 
or set of parameters. In J7] we introduce a new line shape diag- 
nostic and repeat the analysis similar to the EW analysis. Our 
conclusions are summarized in $00 



2. Methodology 

2.1. Stellar Models 

We study the state-of-art atmospher ic models produ c ed by 



the multi-purpose cod e PHOENIX (Hauschil dtet alJ _ 
Hauschildt. Allard. & Barortj, Il999l) . These models are static 
and assume one-dimensional, spherically-symmetric radiative 
transfer. We have chosen PHOENIX because its models can be 
calculated for all stellar parameters from dwarfs to giants, for 
both static and variable stars. A major advantage of PHOENIX 
models is that the code can include self-consistently NLTE ef- 
fects. LTE assumes that, within a confined region, absorption 
and emission of photons occurs at the same rate and the un- 
derlying region's temperature structure follows a black body. 
NLTE models relax this assumption so that atomic physics may 
produce situations in which local areas need not be in ther- 
mal equilibrium. Spectra calculated from the PHOENIX mod- 
els can be calculated for almost any resolution ranging from the 
UV to radio wavelengths. 
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In order to compare stellar model atmospheres to the obser- 
vations of real stars we must "transform" the models into ob- 
servables by simulating the effects of external phenomena on 
the spectra. Microlensing, for example see |5j would be such 
a transformation, as it occurs externally to the source star, but 
does affect the light collected by the observer as a function of 
time. Convolving model data with this transformation produces 
synthetic data, which we can then analyze in the same way as 
we would analyze observational data. We conduct a compar- 
ison between PHOENIX LTE and NLTE synthetic spectra as 
well as a comparison to high resolution observational spectral 
data from a microlensing event. 

3. Ha and its Equivalent Width 

In this work we focus on the effects that can be studied with 
observations of the Ha line. The Ha line is a strong line for 
stars with 4500 < T e ff < 5600/T and is therefore a frequent ob- 
servational target. Different parts of the Ha line are formed in 
different regions of a stellar atmosphere, some which are more 
affected by NLTE than others (Grav, 1992). This makes Ha a 
good candidate to study possible differences between LTE and 
NLTE modeling by testing these differences against observa- 
tions. Another incentive to study Ha is that published data are 
available for t he resolved source star du ring some microlensing 
events (Albrow et al., 2001b; Ca ssan et aill2004l) . In this sec- 
tion, we concentrate on the EW of the line as a figure of merit 
with which to compare LTE with NLTE effects, and models 
with observations. We adopt the usual definition of EW: 



EW 



f 



l - 



Fcont(A) 



dA 



(1) 



where Fn ne is the radiant flux as a function of wavelength A, 
and F cont the value of the flux from the continuous spectrum 
outside the spectral line. The analytical potential of EW as a 
tool for studying atmospheric models depends on one's ability 
to measure its value consistently for both model and observa- 
tional spectra, which requires a careful choice of the contin- 
uum and the rejection of spectral blends. In order to reduce 
systematically-introduced effects to the analysis, we apply the 
same algorithm to all our model and observational data. 

3.1. Extracting Equivalent Width from Model Spectra 
and Data 

For both model and observational data, the EW is determined 
as follows. We first apply a standard vacuum-to-air conversion 
to all model data. We then remove blended spectral lines and 
from this set of model/data points we determine the continuum 
around the spectral line that we wish to analyze and normalize 
the spectra. Note that by "continuum," we mean the "pseudo 
continuum" determined by fitting, as there is no way of deter- 
mining the real continuum for spectral line observations. After 
smoothing the data, we obtain the EW by spline interpolation 
followed by direct integration. 

Our clipping algorithm removes lines blended with Ha 
from the spectral data, which would otherwise contribute to the 
Ha EW during direct integration. We discard all points from the 
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Fig. 1. Illustration of our EW integration algorithm: black solid 
points are the model output data around the Ha line with added 
noise, distributed at random following a Gaussian with a root 
mean squared (rms) scatter <x of 2.7%. Marked in white are 
those points chosen by the algorithm to be used for the inte- 
gration. These are then binned, spline interpolated (black solid 
line), and integrated up to the continuum (long dashed line) to 
give the EW of the line. 

initial spectral data (whether model or fully-reduced observa- 
tional data) that differ in flux from that of either of their neigh- 
bor points by more than a threshold, T, which is determined by 
the signal-to-noise of the data and its spectral resolution, A A 
via, 



T = C<t 2 /AA . 



(2) 



where <x is the rms scatter of all flux points about their mean 
across the wavelength range considered. C is a constant used to 
adjust the magnitude of the threshold to ensure optimal clip- 
ping. We chose C by testing our clipping routine on differ- 
ent modeled and observed spectra, and found that a different 
value of C must be used within the core of the line in order to 
not reject spectral points corresponding to the line itself. After 
clipping, the continuum is then determined from the remaining 
points. 

Figure [0 illustrates how we obtained the Ha EW in both 
synthetic and observed data. We assume the continuum to be 
linear over the narrow spectral range considered. We compute 
the mean flux value of the outermost 1 A of the spectra on ei- 
ther side of the Ha line, and fit a continuum line through these 
values. Although the synthetic data are virtually noiseless, we 
determine the continuum in this way so that the same method 
can be applied to the noisier observational data. We have tested 
and confirmed this method by determining the EW for integra- 
tion ranges from 1A to 10A. The slight increase in EW was 
linear from 9 A onwards, indicating that the line signature lies 
within ± 9 A and it is safe to determine the continuum outside 
of this range. 

When smoothing, we divide the spectral data into three sep- 
arate sections: (a) 6564.8 + 1A to treat the core of the line, (b) 
the region excluding the core out to 6564.8 + 3 A, to cover the 
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wings, and (c) the transition to the continuum of the spectral 
range of 6564.8 + 10A, excluding the core and the wing re- 
gions. Within each of these spectral regions we adapt the sizes 
of bins over which we average so as to sample the continuum 
less densely than the wings, and the wings less densely than 
the core. The averaged flux values in each bin are then spline 
interpolated to give the EW of the line via direct integration. 



4. The Stellar Parameter Grid: LTE vs NLTE 

Using this algorithm, we obtained disk-integrated Ho- EW val- 
ues over a grid of 168 PHOENIX LTE and NLTE models. The 
stellar parameters ranged over a grid of 2.0 < log(g) < 5.0 in 
steps of 0.5 and 4500/T < T eff < 5600/T in steps of 100K. A 
metalicity of [Fe/H]= -0.5 was assumed throughout. 

In Figure|2l we indicate contour lines of equal Ha EW val- 
ues, the results for both LTE and NLTE across the atmosphere 
model grid. 

One can see immediately that Ha EW is very sensitive to 
changes in T^, but much less so for changes in gravity. NLTE 
models with higher r e jf show a slightly stronger dependence 
on gravity than do the respective LTE models. Previous studies 
(e.g. Grav, 1992) have shown that NLTE calculations produce 
a deeper core for the Ha line than in the LTE regime, which 
we confirm. The bottom panel shows the fractional difference 
in Ha EW between the two model types. 

Over the middle range of the grid, the fractional differ- 
ence between LTE and NLTE is about 15%. This shows clearly 
that, within our stellar parameter grid, the difference between 
NLTE and LTE is significant, and that large systematic errors 
may compromise any conclusions drawn for the Ha line in this 
regime if LTE is assumed for its formation. This work also indi- 
cates that over the stellar range T e ff > 54QQK and log(g) < 3.0 
a LTE approach may be reasonable for a study of the Ha line, 
depending on the precision required. 



5. Microlensing as a Test: Resolving The Surface 
Of Stars 

So far we have only compared models of stellar atmospheres. 
Can these differences be observed and, if so, how can this help 
us to understand better the underlying physics of line forma- 
tion and stellar modeling? In particular, do we have ways to 
distinguish between NLTE and LTE effects at the limb and 
the center of stars? We might expect to observe larger differ- 
ences at the limb, where less light originates from deeper lay- 
ers of the star, as NLTE processes are more important in the 
outer, less dense regions of the atmosphere. Resolving the stel- 
lar surface then becomes crucial. A fairly new tool of resolving 
the surface of stars farther away than the Sun is microlensing 
(Brvce, Hendry. & Valls-Gabaud, 2002). When a stellar object, 
the lens, crosses the line of sight to a source star and the lens is 
a binary, an asymmetric magnification pattern, which peaks at 
the so-called caustic, moves across the face of a source star (see 
Figure[3J. During such a caustic crossing there is a strongly dif- 
ferential magnification across the face of the microlensed star, 
which allows its surface to be resolved. 




-2 2 -1 -0.5 0.5 1 

star center to caustic distance [R star ] 

Fig. 3. The effect of a microlensing caustic crossing on the in- 
tegrated light and EW of the Ha line.: Left panel: As a star 
crosses a caustic, it becomes differentially magnified due to the 
asymmetric magnification pattern (red long dashes). Only the 
exit is shown. The vertical thin lines show contours of equal 
magnification leading up to the fold caustic line. Magnification 
drops immediately outside the caustic. Right panels show syn- 
thetic data for which a stellar model atmosphere has been 
convolved with a microlensing magnification pattern. The top 
panel shows how the total flux received from the star changes 
due to the magnification of the caustic. The bottom panel 
shows the respective changes in Halpha EW for this model. 
Note the weaker/stronger dips in EW that occur when the 
leading/trailing limb crosses the caustic, contributing emission 
(negative EW) to integrated light. 

Figure|5]illustrates the microlensing effect on the total in- 
tegrated light of the star (upper right panel) and on the Ha EW 
(lower right panel) during a caustic exit. Light from the parts 
of the star directly coincide with the position of the caustic are 
most highly magnified, and thus contribute a larger fraction to 
the integrated light observed. Figure |4] illustrates the different 
shapes of the Ha spectral line across the face of the star for 
different values of jj — cos 6, where 8 is opening angle to the 
line-of-sight. Ha is a strong absorption line in the center of the 
star (p. — 1) and goes into emission near the limb. This emis- 
sion explains the dip in the Ha EW during the caustic crossing 
while the limb is highly magnified. 

Due to the relative motion of source and lens, the magni- 
fication pattern moves across the source star and produces a 
light curve of integrated stellar light, with a typical shape for 
microlensing events. The magnification pattern may be rather 
complex, but in most cases one can make the assumption that 
the caustic crossing is a fold caustic crossing, ie., that the caus- 
tic has no curvature across the face of the star. For what follows, 
we further assume that the source star is spherically symmet- 
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Fig. 2. Upper panel: contours of equal Ha EW (integrated over the face of the stellar disk) over a range of stellar parameters in 
r e ff and log(g) for LTE atmosphere models. Middle panel: Same, but for the NLTE atmospheric models. For the stellar parameter 
range shown, NLTE models have higher EW compared to LTE models at the same T^ and log(g), as expected for the NLTE- 
affected Ha line. Lower panel: Contours of fractional difference 1 - {EWue/ EWnlte), given in percent. 



ric, and that the relative motion of source star and caustic is 
constant and rectilinear. For fold caustics, a five-parameter fit 
to the light curve data is sufficient to extract the magnification 
pattern of the event and thus infer the contribution each part of 
the source star makes to the integrated light. The five parame- 
ters describe the time and duration of the event, the maximum 
magnification, how long it takes to reach maximum magnifi- 
cation an d the underlying linea r magnification after the caustic 
crossing (Albr owet all 1999b). Entry into a caustic is unpre- 
dictable, but caustic exits can often be predicted so that proper 
observations can be put into place. Outside the caustic, typi- 
cally there is little or no significant differential magnification 
across the stellar disk and the spectrum is similar to that of the 
unmagnified star. During the exit, integrated light is substan- 
tially influenced by the highly magnified light from the limb, 
only where differential magnification is important. 



5.1. Variations of the Ha line EW during a fold caustic 
crossing: A worked example 

Because the Ha line is in emission at the limb, the integrated 
EW of the Ha line drops during a caustic exit compared to 
the non-differentially magnified case. This occurs first when 
the leading limb crosses the caustic and experiences maximum 
magnification. The effect is stronger when the trailing limb 
crosses the caustic, as now most of the stellar disk has crossed 
into the low magnification region outside of the caustic and 
contributes a smaller fraction to the total integrated light than 
when most of the star still lies in the caustic interior. 

During the microlensing event microlensing event OGLE- 
BULGE-2002-069 the PLANET collaboration 2 was able to 
obtai n high re solution VLT UVES data during the caustic exit 
(Cassan et al., 2004), and the collaboration kindly provided us 
with the fully-reduced spectra around the Ha line. 
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Fig. 4. Ha spectral lines from stellar model atmospheres for 
different values of ju - cos 6, where 6 is opening angle to the 
line-of-sight. At the limb // = and fi - 1 in the center of the 
stellar disk. Ha is a strong absorption feature for most of the 
stellar disk, but near the limb Ha goes into emission. 

We adopted t he parameter i zation of the fold caustic exit 
of the event from Cassan et al.l ( 120041) . For each time at which 
a UVES spectrum was available we produced an equiva- 
lent synthetic microlensed spectrum for both LTE and NLTE 
PHOENIX model atmospheres, for each atmospheric model on 
our stellar parameter grid. 

We measured the Ha EW of the synthetic microlensed data 
via the algorithm described in |3] and applied the same algo- 
rithm to the observational data. Due to the virtual non-existence 
of noise in the synthetic data, we adjusted some values of the 
fitting parameters for the observational data: the clipping pa- 
rameter C was increased, and we have added an extra 8A out- 
wardly to the range over which we determine the continuum 
to achieve a higher stability in the continuum. Also, since the 
observed width of the line is not reproduced by the models, as 
we describe later, the "core region" was defined as Ha + 0.9A 
in the data instead of Ha + 0.75A in the models. 

Figure [5] shows our results, where the Ha EW is plotted 
for each time at which observations have been taken. (The 
last two points were taken long after the caustic exit to en- 
sure two baseline points not affected the differential magnifi- 
cation.) We plot the Ha EW measurements for only one set of 
stellar parameters, namely: Tgg - 5000K, log(g) = 2.5 and 
[Fe/H] = -0.5. These are the parameters from our grid of stel- 
lar models th at are closest to the stellar parameters given by 
ICassan et al.l fc(X)4ll J eff = 5000K, log(g) = 2.5 and [Fe/H] 
= -0.6. The difference in metallicity of 0.1 dex is well within 
observational errors of the metallicity determination. Further, 
in previous tests, we did not detect any mesurable changes in 
Ha EW from models with [Fe/H] = -0.5 to models with solar 
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Fig. 5. Changes in the Ha EW for OGLE-BULGE-2002-069 
during the course of its caustic exit. The solid lines represent 
the synthetic data for LTE (upper panel) and NLTE models 
(lower panel) for T eff = 5000K, log(g) = 2.5 and [Fe/H] = - 
0.5, taken fro m our model grid, closest to those published by 
ICassan et alJ ( 120041) . The filled circle s show our d erivation of 
Ha EW from the spectra taken by ICassan et al.l l|2004). We 
corrected the observations for the systematic uncertainties in 
the EW determination for LTE and NLTE models computed 
from Monte-Carlo simulations, and applied their statistical un- 
certainties as errorbars. We did this separately for the high and 
low S/N regions. The double-ringed datapoints show observa- 
tions taken 36 days after caustic exit. 



metallicity, so we are confident in our assumption that there are 
no effects from the 0. 1 dex difference in metallicity. 

We quantify the uncertainty in our determination of the 
Ha EW, using a Monte-Carlo approach simulating 999 noisy 
spectra by adding random noise to the virtually noiseless syn- 
thetic data and rerunning the algorithm to determine the EW. 
The simulations have been restricted to a stellar model of 
r e fj = 5000K and log(g) = 2.5. We assume no variation of the 
fractional error distribution across the whole stellar parameter 
grid. The random noise we add follows a Gaussian distribution 
with an rms scatte r cr equivalent to a signal-to-noise (S/N) of 
130 or 50, which ICassan et al.l d2004h quote for their spectra 
at a resolution of /i/A/l = 30000. These values translate into a 
signal-to-noise per spectral data point, spaced at 0.0174A, of 
(S /N)pp = 36.5, or 2.7% noise, at those times where the cen- 
ter of the star is most highly magnified. At the caustic exit and 
post-exit time we obtain a (S /AOpp = 14, or 7% noise. When 
we add this noise to the model spectra at the appropriate times, 
the model and observed spectra indeed appear to have a similar 
noise level. Thereafter, we shall use high and low S/N or noise 
to refer to the 2.7% and 7% noise levels respectively. 
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Simulation 



^aew[A] o- EW [A] cr EW [%] 



LTE (high S/N) 0.03379 

LTE(lowS/N) 0.03175 

NLTE (high S/N) -0.0623 

NLTE (low S/N) -0.0912 



0.0492 


3.28 


0.1626 


10.8 


0.0586 


3.45 


0.1509 


8.88 



Table 2. Parameters of Gaussian fits to the Monte-Carlo simu- 
lation results for determining systematic and statistical errors. 
We have separately determined these for LTE and NLTE mod- 
els at high and low noise levels Both models have the same 
stellar parameters. The systematic error, //aew is the amount 
that must be subtracted, on average, from the EW determined 
from noisy "data" compared to that determined from noiseless 
models. The scatter in EW determination about its known mean 
iscr EW 



From the 999 simulations per noise level and model 
type, we determined a histogram of AEW = EW mo( jei,noisy - 
EW mo dei,noise-iessfrom noisy to noise-less data, fitted a Gaussian 
to these. The fitting parameters are listed in Table l5.ll 

The scatter in the observational EW points is similar to the 
cr EW determined from the model spectra. Large scatter appears 
to be mainly a continuu m fitting effect . We fi t a continuum to 
the spectra normalized bv lCassan et al . (2004) in the same way 
as we do it for the models. Since the observational data is noisy, 
the continuum is more uncertain compared to noiseless model 
data. Even slight differences in the choice of the continuum 
will result in a larger scatter of E W value s. Since the data have 
already been normalized bv lCassan et alJ (120041) . for the obser- 
vations we fit a continuum line with a fixed slope of zero. Our 
measure ment of the continuum is slightly lower than that cho- 
sen bv lCassan et alJ J2.004f> - but is consistent with the way we 
pick the continuum for the models. 

The EW of the data determined in this way is shown in 
Figure|5]as dots; in the same figure, model predictions are also 
indicated. The overall trend of the observational data is repro- 
duced by the models, however the offsets between them is ob- 
vious. 

The LTE models do not match the magnitude of the EW of 
the Ha line of OGLE-BULGE-2002-069 after the caustic exit 
(double-ringed data points), where no differential magnifica- 
tion occurs and also at times when the center of the source star 
crosses the caustic. This could indicate that either the stellar 
atmosphere models are insufficient and/or we have chosen the 
wrong stellar parameters from the grid. 

6. Giving the models a fair chance 

Before condemning the model atmospheres prematurely, we 
need to be sure that we have chosen models that best fit the 
Ha EW data within uncertainties in the determination of T e g 
and log(g). We quantify the discrepancy between models and 
observations by defining^- 2 , the goodness-of-fit figure of merit. 

X 1 = J] [EW obs (tobs) - £W mode i(f obs )] 2 /(<r 2 w (f obs ) , (3) 

fobs 

summing over all points in time f b s at which an observed spec- 
trum has been taken. We assume that any systematic errors in- 



troduced by our fitting algorithm and statistical uncertainties in 
fitting model spectra are the same across the model grid. Thus, 
for ctew we use the values derived in 35.11 namely, around 10% 
and 3.5% for the high and low noise data, respectively. 

We calculate x 1 f° r eacn model in our model stellar grid, 
and show in Figure|6] contours of equal Ay 2 , 



A* 2 =* 2 



■x 2 - 

A mm 



(4) 



where v 2 . is the minimal value of v 1 across the model grid, 
i.e. the goodness-of-fit for the best fitting model. Models along 
a contour line fit the observed measurements equally well. In 
Figure[6] we superimpose on the model grid the stellar param- 
eters determined by ICassan et al. (2004 ) for OGLE-BU LGE- 
20 02-069. We follow a n analysis of lMinnitietaD(l2002i) . used 
bv lCassan et alJ ( l20 04). indicating uncertainties in determining 
stellar parameters of +200K in T e ff, +0.5 in log(g), and ±0.3 in 
[Fe/H]. 

Within the uncertainties of the stellar model, LTE models 
with x 10QK hig her 7W than and NLTE models at the same 
values quoted bv lCassan et alJ (120041) provide the smallest^ 2 . 
Gravity is not well constrained. The overall goodness of fit is 
similarly poor for both LTE and NLTE models at a values of 
« 28 for 20 observations, which could be caused by random 
noise, with only about 7% confidence. Only the best fitting 
models should be examined to study further differences to the 
observations. Note that the best-fitting LTE and NLTE models 
have different temperatures. 

7. The Importance of Line Shape 

Equivalent width collapses all the information about a spectral 
line into one value that only considers the area of the line under 
the continuum, but is independent of the line shape. This means 
that, despite a good agreement between the EW of a model 
spectral line, and of the observed spectral line any variations in 
the line shape will be ignored if this is the only figure of merit 
used to test models. 

Examples of the Ha spectral line (top panel) and their fit 
from the algorithm (bottom) panel are shown in Figure [7] for 
one observation. The line from corresponding differentially 
magnified LTE and NLTE models with parameters as chosen 
from ICassan et alJ (2004) are shown in the left panels, right 
panels give the best fitting stellar models as determined in 
Figure |6] Although the changes appear small, allowing stellar 
parameters T e ff and log(g) to vary within the observationally al- 
lowed uncertainties, can alter the EW by 10-15%, as can been 
seen in Figure|3 

That the models do not generally fit the width of the Ha 
line (see Figure [7) suggests that there are other physical pro- 
cesses occurring which are yet to be explained by stellar mod- 
els. One e xplanation may be chromospheric influences, which 
ICassan et alJ ([2004) claim to have detected in their data. Since 
Ha is a chromospheric line but the standard model atmospheres 
do not include chromospheric calculations, more research is 
needed to explain this phenomenon. 

This comparison of individual spectra clearly shows where 
discrepancies in EW occur, primarily in the wings of the line, 
but for LTE models also in the core of the line. However, the 
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Fig. 6. Upper panel: Contours of equal Ay 2 - X 2 ~ X 1 ■ > where Ay 2 . is minimal when best agreement between LTE models 
and observations of the Ha EW of OGLE-BULGE-2002-069 are achieved. Lower panel: Same for NLTE. A% 2 is different for 
LTE and NLTE models and its value is shown in the left corner of the figure panels. The grey area indicates all models that lie 
within 3<x of the best fitting model for two degrees of freedom, namely T e g and log(g). Further contour lines are the 4<x, 5<x and 
highe r boundaries. The tick marks on the contour lines indicate the direction of decrease in A;^ 2 . From their spectral analysis, 
ICassan et alJ ( 120 04) determined the stellar parameters to be T e g - 5000K, log(g) = 2.0, [Fe/Hl= -0.6. We mark this c hoice of 
parameters with the large cross with uncertainties of AT e g - +200K and log(g) = +0.5 as given in Minniti et alJ ( 120021) for this 
type analysis. 



differences in shape in the core and the wings balance each 
other out when only total EW is considered. The discrepancies 
due to the different line shapes are washed out. Therefore we 
have devised a different diagnostic to compare model spectra 
with observational data by defining the line shape diagnostic, 



(A (fobs) 



«71ine 



[^obs,o(^, fobs) - 7model,o(4, fobs)] dX 



where we integrate the square of the flux difference that re- 
sults when we subtract the model spectrum from the observa- 
tional spectrum, after both have been normalized, clipped and 
spline-interpolated. Note, that the units of this new diagnostic 
are A, since the flux has already been normalized. Similarly, 
we can obtain if/ 2 for line shape comparisons between LTE 
and NLTE models. The result of which is shown in Figure [S] 
for non-microlensed model data. In order to compare LTE and 
NLTE models to observational data taken during a microlens- 
ing event, and to find the best-fit model based on the measured 
shape parameter, we define a new shape figure of merit 




TsJ]^«oh) 



(6) 



which is sum of equally weighted iff 2 over each observation 
during the event. The units of *P are still A. 



45 46 47 48 49 50 51 52 53 54 55 56 
T.„ [100K] 

Fig. 8. Contours of equal tfr 2 for a NLTE-LTE model compar- 
ison. The tick marks on the contour lines indicate the direc- 
tion of decrease. Measured by this quantity, line shapes for 
NLTE and LTE models differ least for lower r e ff and higher 
log(g) within the parameter grid. Larger differences are found 
for higher T e ff and higher log(g). 



Defining the shape parameter *F allows us to collapse the 
available time series data to one quantity per stellar model, 
while conserving the information of the line shape compar- 
ison. In effect, Y is sum of x 1 comparison of all the spec- 
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Fig. 7. Differences in the spectra between observations of the Ha line for OGLE-BULGE-2002-069 (in black) and models (LTE 
in blue, NLTE in red) at the same magnification. The left panel shows comparison to the models with the stellar parameters 
closest to the published values, the right panel shows the spectra from models as chosen from Figure [6] Bottom panels present 
the clipped and spline interpolated lines after the observations and models have passed through the EW algorithm from 33.11 
Each inset shows a narrower region around the core of the line. For all diagrams log(g) = 2.5 and T e s is given per 100K. 



tra themselves, across a wavelength range defined by the line, 
equally weighted. We determine uncertainties in *P by fitting 
a Gaussian to the distribution of *P values resulting from 999 
Monte-Carlo simulations. In each simulation we sum up 13 
low noise and 7 high noise if/ 2 values randomly drawn from 
the distribution of if/ 2 values, which we have determined via 
999 Monte-Carlo simulations in a similar fashion as for EW in 
35.1l for low and high noise, LTE and NLTE models. We fit the 



resulting distributions of if/ 2 with a combination of a Gaussian 
and a straight line as the distribution is asymmetric because if/ 2 
is positive definite. We fit the *T distribution with a Gaussian 
convolved with a straight line. From this we deduce 3, 2 and 
lcr levels of confidence which we mark in Figure [9] as dark 
grey, light grey and the enclosed white areas. 

In Figure|5] we show contour plots of equal T'-TVin, where 
Tmin is the lowest value of ¥ across the model grid. For LTE 
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Fig. 9. Contours of equal *F - *P m i n for comparison of the observational data with LTE models (upper panel) and NLTE models 
(lower panel). The contours clearly constrain not only the T e ff but also the log(g) of the best fitting model parameters. This means 
that the best fitting LTE model has not o nly a higher 7W a nd a lower log(g), whereas the best fitting NLTE model has stellar 
parameter closest to the one suggested bv lCassan et alJ l|2004), indicated by the black cross. 



models this is ^ min = 0.71, for NLTE models ¥ min = 0.47. The 
shape parameter if/ 2 is much more dependent on log(g) than is 
EW, which was strongly dependent only on T e g. The dark grey 
shaded area indicates the 99.73% level of confidence around 
the best fitting model. The enclosed light grey and white ar- 
eas show the 95.4% and 68.3% confidence levels, respectively. 
The best-fitting LTE model has a log(g) of less than 2.0, sub- 
stantially less than what the observationally-determined stellar 
parameters suggest. As with fits to EW, LTE fits to the shape 
parameter, V F, again suggest a higher r e ff by about 200K, i.e. 
100K higher than the EW-best-fitting LTE model suggested. 
Although the best-fitting NLTE model for the shape parameter 
leans towards lower gravities, these are in agreement with the 
observationally determined stellar parameters. The overall fit 
to *P, however, is fairly poor: The distribution of *¥ at random 
shows that 99.9% of values drawn randomly from the distribu- 
tion of if/ 2 for NLTE models suggests a value below 0.07. For 
LTE models this reduces to a value less than 0.01. 



8. Conclusions 

As a part of a larger program to compare model atmospheres 
to observational data, we have developed an algorithm to con- 
sistently and automatically determine and compare spectral 
line characteristics. Here, we have focused on equivalent width 
(EW) of the Ha line and also introduced a line shape parameter 
if/ 2 , which emphasizes flux differences between two spectra. We 
have studied a stellar parameter range of 4500 < T e g < 5300 



and 2.0 < log(g) < 5.0 at [Fe/H]= -0.5, a range typical of 
microlensed source stars in the Galactic Bulge. 

We have determined the uncertainty of this algorithm by 
performing Monte-Carlo simulations in which noise is intro- 
duced to a model spectrum and the extraction procedure re- 
peated. The simulations have been restricted to a stellar model 
of T e s = 5000K and log(g) = 2.5. We assume no variation 
of the fractional error distribution across the whole stellar pa- 
rameter grid. The EW of the Ha spectral line as predicted by 
PHOENIX atmospheric models is substantially different for 
LTE and NLTE model. The wrong choice of model type may 
result in uncertainties in EW larger than 8% for most models 
with 4500 < r eff < 5300 and 2.0 < log(g) < 5.0. The shape 
parameter if/ 2 in an LTE-NLTE comparison is statistically sig- 
nificant and largest for high T e ff and high log(g). 

Choosing the right model to compare to the observations 
is crucial when one wants to understand stellar atmosphere 
modeling. As an example, we study Ha spectral data from the 
microlensing event OGLE-BULGE-2002-069. First, we deter- 
mine the best-fitting model to the Ha EW data over the range of 
our stellar grid for both LTE and NLTE models over the course 
of the caustic exit, during which the source star is resolved. We 
find that in order to fit the observed data an LTE model requires 
a r p ff about 100K larger than that proposed by ICassan et alJ 
(2004) based on spectral type analysis for this source star. The 
best-fitting N LTE model, on the other hand, requires a T e ff as 
proposed bv lCassan et alJ ( 120041) . The absolute x 1 is lower for 
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NLTE PHOENIX models than for LTE models, suggesting a 
slight preference for the NLTE model. 

Second, we determine the best-fitting model to the Ha if/ 2 
which we define, and find that the best-fitting LTE model has a 
J e ff, that is 200K higher than the observationally determined 
r e ff and a significantly lower log(g), The best-fitting NLTE 
model agrees with the observationally determined values for 
r e ff and log(g) within the uncertainties, although it tends to un- 
derestimate both T e ff and log(g) for this star. The NLTE models 
produced a better figure of merit than the LTE models for both 
diagnostics. 

We caution, however, that declaring either LTE or NLTE 
models a failure based on analysis of Ha only is premature in 
the case of OGLE-BULGE-2002-069 . Chromospheric effects 
may alter the characteristics of the Ha line significantly, and the 
atmosp here models used h ere do not include chromospheres. 
Indeed. lCassan et alJ (120041) claim to have detected an anomaly 
in the Ha line at some times during the caustic crossing that 
may be due to a chromosphere. 

We would like to point out that with our method we are 
already placing strong constraints on T e ff and log(g) by using 
only two parameters measured for one particular line. Future 
studies of multiple lines and multiple diagnostics will provide 
even stronger constraints. In subsequent work, we will deter- 
mine the most effective choices of spectral lines and their di- 
agnostics for testing model atmospheres and include chromo- 
spheric calculations in the models for direct comparison to ob- 
servations of OGLE-BULGE-2002-069. 
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